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Abstract. 

By a few but important examples as models of combined radiative and 
wind-driven Hii regions and galactic winds we demonstrate the importance of 
refined small-to-medium scale studies of chemo-dynamical effects. These pro- 
cesses determine the internal dynamics and energetic of the ISM and affect its 
observational signatures, e.g. by abundance contributions, but are not yet reli- 
ably and satisfactorily explored. 



1. Introduction 

Galaxies experience energy deposit and element release into their ISM from 
massive stars during the lives and by their explosions. They are main drivers of 
energetics and dynamics in galaxies and in both stages they release characteristic 
elements which are preferably contained in the hot gas phase. Nevertheless, 
because of the close coupling of hot, warm, and cool ISM energetic and dynamical 
processes lead to chemo-dynamical effects which alter the energy budget and the 
element content in the different gas phases. Understanding these processes from 
stellar to galactic scales are therefore vital for understanding and modeling the 
evolution of galaxies and also their interaction with the environment. 



2. Massive Stars 

2.1. Effects of massive stars during their lives 

Massive stars play a crucial role in the evolution of galaxies, as they are the 
primary source of metals, and they dominate the turbulent energy input into 
the interstellar medium (ISM) by their massive and fast stellar winds, the ultra- 
violet radiation, and supernova explosions. The radiation field of these stars, at 
first, photo-dissociates the ambient molecular gas and forms a so-called photo- 
dissociation region (PDR) of neutral hydrogen. Subsequently, the Lyman con- 
tinuum photons of the star ionize the Hi gas and produce a Hii region that 
expands into the neutral ambient medium. 

As these stars have short lifetimes of only a few million years, Hii regions 
indicate the sites of star formation (SF) and are targets to measure the current 
SF rate in a galaxy. Furthermore, the emission line spectrum produced by the 
ionized gas allows the accurate determination of the current chemical composi- 
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tion of the gas in a galaxy. Although the physical processes of the line excitation 
are quite well understood and accurate atomic d ata are available, so that the 
spectral analysis of Hii regions (see e.g. [Stasiiis ka 19791 : lEvans Sz Dopita 19871 ) 
serves as an essential tool to study the evolution of galaxies, their reliability as 
diagnostic tool hav e also to be studied with pa rticular emphasis e.g. to temper- 
ature fluctuations ( Peimbert 19671 : IStasihska 2002 ) and line excitations. 

While the simple concept of a uniform medium in ionization equilibrium 
with the radiation from the stars (the Stromgren sphere) is successful in de- 
scribing several global features of the object and allows to predict the emission 
line spectrum quite reliably, it has long been realized that Hii regions are com- 
plex and dynamical objects. Dynamical modelling of Hii regions caused purely 
by the energy deposit of the stellar radiation field has therefore been started 
already long ago (see e.g. lYorke 1986I . and references therein) providing both a 
first insight into the formation of dynamical structures and allowing to derive 
more realistic scaling laws for observational comparisons, respectively. 

That massive stars do not only release energy by their radiation but also act 
by means of strong stellar winds which drive a turbulent ISM but, on the other 
hand, both are metal dependent has three major consequences for galaxies with 
present-day solar metallicities: at first, the classical Hii has lost its simplicity, 
because, secondly, it engulfs a vehemently expanding stellar wind bubble (SWB) 
and, thirdly, peels off the stellar outermost shells, by this, incorporating nuclear 
burning regions into the wind mass loss. 



2.2. Star formation Feedback 

Quiescent SF as a self-regulated process is a widely accepted concept. The stel- 
lar feedback can adapt both signs, positive as a triggering mechan i sm in a self- 

propagating manner like in superbubble shells flEhlerova et al. 19971 : iFukuda &: Hanawa 2000l ) 



vs. negative as self-regulation. Primarily the correlation between the sur- 
face density of disk gal axies' Hi gas and the vertically integrated SF rate de- 
rived from the Ha flux ( Kennicutt 19981 ) serves as the best proof of a SF self- 



regu lation. Even more refi ned observations that trace more complex molecular 
gas ( Gao &: Solomon 20041 ) follow the same surface density relation. Surpris- 



ingly, this relation can be reproduced by different numer ical models an d with dif- 
ferent prescriptions from pur ely collapsing cloud cores (|Li et al. 20051 1 via ther- 



modynamical cloud models (jKrumholz &: McKee 20051) to sophisticated global 



chemo-dynamical treatments of galaxy models ( Burkert et al. 19921 : iHarfst et al. 200^ ) 
and even seems to hol d under the neg lection of self-regulation processes in cos- 
mological simulations ( Kravtsov 20031 ). This fact obtrudes that the gas depen- 



dence of processes within the SF matter cycle although expected to behave 
non-linearily cancels out and results in a universal law. Such a dependence is 
e.g. plausible because Koppen et al. (1995) have demonstrated that due to the 
square-dependence of collisionally excited radiative cooling on the gas volume 
density also the SF rate adapts to that over a large density range and almost in- 
dependent of the heating rate and other flducial parameters. This also holds for 
a multi-phase ISM accounting for supernova typell (SNII) p roduction of hot gas 
and for gas transitions by evaporation and condensation (K oppen et al. 19 98). 
From this result it seems reasonable that the vertical accumulation of gas mass 
should also lead to a SF rate dependence on the column gas density because the 
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self-gravitational potential determines the disk stratification. The Kennicutt law 
therefore could reflect a sum over the volumes, while the 1.4 power is also close 
to the pure free-fall criterion.. Then the timescales must definitely tell us about 
both process and carrier of self-regulation, respectively. 



2.3. Evolution and structure of bubble around massive stars 

However, this question needs more exploitation yet. Already the energy deposit 
by massive stars described by a parameter e, called energy transfer efficiency, 
is neither well derived during their normal lives nor for their SNII explosions 
and the expanding SN remnants. The reason is that for both stages the envi- 
ronmental state and thus preceding processes as well as dynamical effects play 
a significant role. Analytical est imates for p urely radiative Hii regions yield 
e of the order of a few percent ( Lasker 19671 ). Although the additional stel- 



lar wind power L,„ can b e easily evaluated from model and observational values 
( Kudritzki Sz Puis 20001 ). its fraction that is transferred e.g. into thermal energy. 



i.e. the thermal e, or into turbulent energy is not obvious from first principles. 
The picture is that the fast stellar wind creates shocks that form the SWB filled 
with very hot plasma, which expands into the Hii region so that also this Hii 
shell has to expand into the surrounding ISM due to its overpressure. 

The structure and evolution of SWBs and the t ransition to its surr ound- 
ing Hii region can be described by a set of equations ( Weaver et al. 19771 ) from 



which the wind-energy deposit can be derived theoretically under the simpli- 
fying assumptions of a point source of a constant and spherically symmetric 
strong wind that interacts with a homogeneous am bient ISM. The e's amount 
to significant fractions of (see Frever et al. 20031 . and references therein) . It 



has become evident that the stellar parameters such as the mass-loss rate, the 
terminal velocity, the effective temperature and the luminosity of the star vary 
strongly during the stellar lifetimes. While most previous studies dealt with 
either the evolution of Hii regions or of SWBs separately and with constant 
values, little is known about the interaction of these two structures. 

Although the analytical and semi-analytical solutions for the evolution of 
SWBs have been improved over the years as well as the numerical simulations 
have be en done with inc reasing complexity, like e.g. 2D calculations of SWBs (see 
e.g. bv lRozvczka 19851 . and a series of papers) and/or combined ID radiation- 
hydrodynamical models of Hii re gion coupled with the dynamical SWB (e.g. by 



Garcia-Segura Sz Mac Low 19951 . for references see Freyer et al. 2003), a variety 
of physical effects remains to be included in order to achieve a better agreement 
of model s and observa t ions with regard to the evolution of the hot phase in 



bubbles (|MacLow 20001 : IChu 20001 ^ 



To improve the insight into the evolution of radiation-driven + wind-blown 
bubbles around massive stars, we have performed a series of radiation- hydrodynamics 
simulations with a 2D cylindrical-symmetric nest ed-grid scheme for stars of 
masses 15 MfD (Kr oeger et al. i n prep.), 35 Mp, ([Frever et al. 200^ ). 60 M© 
(|Frever et al. 20031 ). and 85 Mq (iKroeeer et al. 20071 1. The main issues are: 



1) The Hii region formed around the SWB has a complex structure mainly 
affected by dynamical processes like e.g. shell instabilities, vortices, mixing 
effects, etc. It is more compresssed by stronger winds (see Fig.l). 
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Figure 1. Temperature distributions of a 15 M© (left panel) and a 85 Mq 
(right panel) model at late stages of their stellar evolution. The temperatures 
are colour coded to the right of the plots but differently (taken from Kroeger 
2006). For discussion see text. 

2) Finger-like and spiky structures of different densit ies and temperatures are 
formed in the photo-ionized region (|Frever et al. 20031 ). 



3) The regions contributing to the Hll emission line are not solely limited to the 
photo-ionized shell around the SWB but also form from photo-evaporated gas 
at the trailing surface of the SWB shock front (see 85 Mq model in Fig.l). 

4) Because dispersion of this cooler photo-evaporated gas into the hot SWB 
leads to mixing also the stellar material expelled by the wind emerge partly in 
the spectra. 

5) As a consequence the metal-enrichment of the wind in the Wolf-Rayet stage 
which is generally assumed to remain only in the hot SWB for a long time 
affects the observatio nally discernible abundance of the Hii gas. By these models 
( Kroeger et al. 20061 ) it could be ])roven for the first time that the metal release 



by Wolf-Rayet stars can be mixed within short timescales from the hot SWB 
into the warm ionized gas and should become observationally accessible. As the 
extreme case for the 85 Mq star we found a 22% enhancement of Carbon, but 
neglible amouts for N and O. 

6) As expected from the distribution of Hii gas the radially projected Ha bright- 
ness shows a decrease to the center and a slight brightening to the limb but not 
as strong as ex pected according t o the increase of the line-of-sight with im- 
pact parameter ( Frever et al. 200^ ). This effect depends on the bubble age and 



starts from central brightning. It also demonstrates both: the neglection of heat 
conduction and the homogeneous initial density do not allow a sufficient bright- 
ning of heat conductive interfaces so that, secondly, only the photo-evaporated 
backflow can contribute to the Ha luminosity in present models. In reality, 
condensations which become embedded into the hot SWB are exposed to heat 
conduction. 

7) The sweep-up of the slow red supergiant wind by the fast Wolf-Rayet wind 
produces remarkable morphological structures and emission signatures which 
agree well with observed X-ray luminosity and temperature as well as with the 
limb brighteni ng of the radially projected X-ray intensity profile (for details see 
iFrever et al~0 06). 

8) e's for both radiative as well as kinetic energies remain much lower than 
analytically derived (more than one order of magnitude) and amount to only a 
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few per mil ( Hensler 20071 ). There is almost no dependence on the stellar mass 
what is principly expected because the energy impact by Lyman continuum 
photons and by wind luminosity increase with stellar mass. Vice versa, since 
the gas compression is stronger by a more energetic wind also the energy loss 
by radiation is more efficient. 

Nonetheless, a word of caution and unfortunately of discouragement has to 
be expressed here because the stellar evolutionary models are not unique but 
depend on the authors. In order to get a quantifiable comparison of the models 
by Garcia-Segura et al. (1996a, 1996b) with our 2D radiation-hydrodynamical 
simulations (Preyer et al. 2003, 2006) for the 35 and 60 Mq studies we used 
the same stellar parameters. Since no stellar parameters were available from the 
same gr oup for the 15 and 85 Mq models we had to make use of the Geneva 
models (|Schaller et al. 19921 1. A comparison of the age-dependent 60 Mq pa- 
rameters of both groups has revealed enormous differences in the energetics by 
almost one order of magnitude as well as that the Wolf-Ray et and Lumiii ous 
Blue Variable stages occur in contrary sequence, respectively ( Kroeger 20061 ). 

Furthermore, Hii regions show a great diversity of shapes which give direct 
evidence that the ISM is not a smooth, homogeneous, and uniform medium. In 
addition, the complexity of the internal structure and of the flow pattern of the 
gas indicates that the evolution of the ionized region is strongly affected by exist- 
ing irregularities and hydro dynamical instabilities. And even more, stellar wind 
bubbles and SNRs of neighbouring stars collide and may overlap to form super- 
bubbles with diameters of several 100 pc scale, which are prominent by their X- 
ray emission in galaxies with active SF. This picture of SWB/Hii aggregates has 
also to be extended by PDRs which are mainly produced in the early evolution- 
ary phases of ultra-compact Hii regions. Si nce they are heavily du st-enshrowded 
they emerge from GLIMPSE observations ( Benjam in et a l. 20031) and allow de- 
tailed diagnostics of these common stellar bubbles (|Chruchwell et al. 20061 ). 



3. Galactic Mass Exchange 
3.1. Galactic Outflows 

Finally, massive stars explode as SNell, creating a SN remnant that sweeps up 
the ambient medium. These phases have been understood and modeled in de- 
tail energeticly and dynamically. Cumulative SNell form superbubbles which 
expand vehemently and can drive a galactic outflow. Since the hot superbub- 
ble gas carries SNII elements and has to act dynamically with the surround- 
ing ISM, it has huge energetic and chemial effects on galaxies. Therefore, the 
questions that arise for their chemo-dynamical contribution to galaxy evolu- 
tion are the same as before for th e Hil /SWB regions. Although investigations 
have been performed for e of SNe ( Thornton et al. 19981 ) and superbubbles (e.g. 



Strickland et al. 20041 ) they are yet too simplistic (e.g. only ID) for quantitative 



results. 

Because of their low gravitational potential dwarf galaxies are sensitive 
to energetic processes like superbubbles so th at these can easily d evelop to a 



galactic outflow as p erceivible e. g. in NGC 170 5 (jHensler et al. 19981 ). NGC 1569 



([Martin et al. 20021 ). I Zw 18 (|Martin 199 "6|) , and many other Blue Compact 



Dwarf Galaxies (BCDs). Dynamical models (jMacLow &i Ferrara 19991 ) can still 
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provide only a qualitative understanding. Nevertheless, the element abundances 
stemming from different stellar progenitors should enable one to trace back the 
history of element enhancement vs. depletion by different processes. Since in 
general the SF varies in BCDs and, furthermore, produces massive star clusters 
one can imagine that a first superbubble has to exert most of its energy for 
its expansion while the successive SNII explosions can be funneled through the 
chimney and by this more easily expel its metal-enriched gas. This can only 
happen, if the next SF episode follows shortly enough after the former one. For 
this scenario it is, thus, important to explore the timescale of refill superbubbles 
vs. the S F timescale. For the refi ll of superbubbles in the gas layer of typical 
BCDs we ( Recchi Hensler 20061 ) found a range of a few 100 Myrs. This means, 
that the amount of galactic winds depends not only on the SF rate but also on 
i ts modes, bursting vs. p asping, but affects the chemo-dynamical evolution 
(iHensler fc Rieschick 20a 
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Figure 2. Comparison of abundance evolution in 2 modls representing 
NGC 1569: NGC-3 with a single strong and short SF burst, NGC- 
4 with 2 former SF epochs and a mild recent burst. (For details see 
iRecchi fc Hensler 2006( ). 



For the best s tudied BCDs, I Zw 18 (|Recchi et al. 20041 ) and NGC 1569 
( Recchi et al. 20061 ). 2D dynamical models with characteristic tracer elements 
have been performed aiming at achieving a best fit with observations. The 
main results can be briefiy summarized as follows: Galactic winds always occur. 
Models with bursting SF are generally unable to account for the chemical and 
morphological properties of NGC 1569, because they either underproduce O or 
inject too much energy into the system, enough to loose a too large gas fraction. 
The best agreement with the chemical composition is found for long-lasting con- 
tinuous episodes of SF of some hundred Myrs of age and a recent more intensive 
short SF burst (see Fig. 2). In most models with gasping SF, the final chemical 
composition of the galaxy reflects mostly the chemical enrichment from old stel- 
lar populations. In fact, if the first episodes of SF are powerful enough to create a 
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galactic wind or to heat up a large fraction of the gas surrounding the SF region, 
the metals produced by the last burst of SF are released into a too hot medium 
or are directly ejected from the galaxy through the wind. They do not have the 
chance to pollute the surrounding medium and to contribute to the chemical 
enrichment of the galaxy. This result confirms the concl usion that most of the 
O from the last SF epoch is stored in the hot X-ray gas ( Martin et al. 20021 ). 



3.2. Gas Infall 

Although the element abundance properties of most BCDs favour only a young 
stellar population of at most 1-2 Gyrs, they mostly consist of an underlying old 
population as it is the case for NGC 1569. Furthermore, most of the objects are 
embed ded into Hi envelops from which at le ast NGC 1569 definitely suffers gas 
infah (|Stil fc Israel 20021 : iMiihle et al. 20051 ). This has lead Koppen & Hensler 



(2005) to exploit the influence of gas infall with metal-poor gas into an old 
galaxy with continuous SF on particular abundance patterns. Their models 
could match not only the observational regime of BCDs in the [12-|-log(0/H) 
-log(N/0) ] space but also explain the shark-fln shape of their data distribution 
by means of the ratio of mass infall to existing cloud mass in the sense that for 
larger infall fractions the loop is more extended. 

Analytical investigations as well as too simplified numerical models, how- 
ever, are unable to distinguish between the different gas phases and their metal 
content in order to derive abundances in the targets of observations. 

In addition, also mixing effects at interfaces between gas phases and due to 
turbulence (like e.g. those in the combined SWB/Hii complexes) contribute to 
the observations by the enhancement or, respectively, dilut ion of metals. Since 
gas infall is not only affecting the chemistry but also SF ( Hensler et al. 20041 ) 



and dynamics of outflows we h ave extended the forme r models (Recchi et al. 
2004, 2006) by infalling clouds ('Recchi fc Hensler 20071 ). In these models, two 
kinds of cloud contributions are considered, only initially existing and contin- 
uously formed, respectively. The issues are the following: Due to dynamical 
processes and thermal evaporation, the clouds survive only a few tens of Myr. 
The internal energy of cloudy models is typically reduced by 20 - 40 % compared 
to models with a smooth density distribution. The clouds delay the develop- 
ment of large-scale outflows, helping therefore in retaining a larger amount of 
gas inside the galaxy. However, especially in models with continuous creation 
of infalling clouds, their ram pressure pierce the expanding supershells so that 
through these holes freshly produced metals can more easily escape and vent 
into the galactic wind. Moreover, assuming for the clouds a pristine chemical 
composition, their interaction with the superbubble dilute the hot gas, reducing 
its metal content. The resulting final metallicity is therefore, in general, smaller 
(by ~ 0.2 - 0.4 dex) than the one attained by smooth models. 
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